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IV 


1.  Introduction 


The  theory  of  phase  equilibrium  for  the  “liquid/vapor”  phase  transformation  has  been  under 
development  for  centuries.  In  the  case  of  one-component  phases,  it  is  presented  in  a  very  clear 
form  in  multiple  classical  monographs  and  textbooks. 1-3  The  key  products  of  the  classical  theory 
of  liquid/vapor  phase  co-existence  are  chemical  potentials  p  of  the  phases  and  the  phase 
equilibrium  curve  p  =  p  (T).  The  phase  equilibrium  curve  gives  the  pressure  p  at  equilibrium  as  a 
function  of  the  absolute  temperature  T.  This  curve  is  universal  in  the  sense  that  it  does  not  depend 
upon  the  mechanical  and  thermal  loading  conditions.  In  other  words,  it  is  one  and  the  same  for 
isothermal  and  adiabatic  systems  as  well  as  for  isochoric  or  isobaric  loading  conditions.  In  short, 
the  equilibrium  pressure  p  depends  upon  the  absolute  temperature  T  alone.  In  the  work  that 
follows,  we  consider  the  liquid  and  vapor  phases  as  co-existing  without  co-mingling.  That  is,  the 
two  phases  are  in  contact  and  may  transport  mass  and  heat  across  a  notional  divide  that  separates 
them. 

Consider  a  substance  with  the  specific  internal  energy  density  e  =  e  (p,  rf),  where  p  is  the  mass 
density  and  77  is  the  specific  entropy  density.  In  order  to  find  the  function  p  (T),  one  has  to  use  the 
following  three  equilibrium  conditions: 

Pi  =  P2,  Ti  =  T-2,  Pl=  p,2  ,  (1) 


where  1  and  2  refer  to  the  liquid  and  vapor  phases,  respectively. 

The  first  of  the  relationships  in  Eq.  1  reflects  the  mechanical  equilibrium  between  phases,  the 
second  reflects  the  thermal  equilibrium  between  the  phases,  and  the  third  reflects  the  “chemical” 
equilibrium.  In  this  context,  the  chemical  equilibrium  means  the  equilibrium  with  respect  to  mass 
exchange  between  the  phase  in  the  process  of  vaporization/condensation. 


Thermodynamical  identities  allow  one  to  express  p,  T,  and  //  in  terms  of  the  derivatives  of  the 
internal  energy  density  function  e  (p,  rf): 


2  de{p,rf)  2  ,  \  _  de  (p,  rf) 

p  (p,  v)  =  p  — Yp —  =  p  ep  T  w  v )  =  — drj — 

P  m 

p  =  e-\ - f-  Trj  —  e  +  pep  +  rjev 

P 


Analysis  of  the  system  1,  which  can  be  found  in  any  thermodynamics  textbook,  results  in  the 


1 


following  equation,  known  as  the  Clausius-Clapeyron  equation:1- 2 

M T)  Vi  -  V2  L 

dT  pi 1  -  p2 1  t  (pr1  -  P2 x) 

The  quantity  L  =  T  (7/1  —  r/2)  is  called  the  latent  heat  of  phase  transformation  and  appears  to  be  a 
state  function,  not  a  process-dependent  function.  In  other  words,  the  latent  heat  L  depends  only 
on  the  thermodynamic  state,  regardless  of  whether  this  state  was  achieved  by  way  of  isochoric, 
isothermal,  or  any  other  particular  process. 


To  make  our  point  regarding  L  more  clear,  let  us  compare  it  with  the  similar  concept  of  the 
specific  heat  C.  Consider  any  process,  in  which  the  heat  supply  0  and  the  absolute  temperature  T 
are  functions  of  time  t  :  Q  —  Q  (t) ,  T  —  T(t).  Then,  the  heat  capacity  C  of  this  process  is  a 
function  of  time  also  C  —  C  (t).  Here,  the  heat  capacity  C  (t)  in  the  process  is  defined  as 


=  dQ(t)/dt 
{  dT(t)/dt 


(4) 


Obviously,  the  function  C  (t  )  essentially  depends  on  the  process  under  consideration.  During  an 
adiabatic  process  it  is  equal  to  zero,  whereas  for  an  isothermal  process  it  is  equal  to  infinity,  etc. 
Using  the  First  and  Second  Laws  of  Thermodynamics,  we  get,  in  addition  to  Eq.  4,  the  following 
formula  for  the  heat  capacity  in  an  isochoric  process,  in  terms  of  the  specific  free  energy  function 

4>(P,T)  =  e-pT: 


Cv 


rdMp,T) 
dT 2 


=  T 


(5) 


Although,  originally,  the  function  Cv  was  introduced  by  formula  4  applied  to  the  isochoric 
process,  in  fact,  according  to  formula  5,  it  is  a  function  of  thermodynamic  state  (p,  T)  only.  In 
other  words,  it  can  be  introduced  directly  with  formula  5  and  without  mentioning  of  any 
particular  thermodynamic  process. 


In  addition  to  Cv  one  can  introduce  the  specific  heat  Cv  for  the  isobaric  or  any  other  specified 
thermodynamic  process.  That  would  lead  us  to  other  analogies  of  Eq.  5  and  additional  state 
functions. 


If  we  introduce  the  state  function  Cv  directly  with  the  help  of  identity  5,  we  cannot  help  asking 
the  question,  “In  what  thermodynamic  process  will  the  function  T  iptt  be  equal  to  the  heat 
capacity  of  the  process?”  In  our  case,  the  answer  is  almost  a  trivial  tautology:  “The  function 
T  ipTT  gives  the  heat  capacity  in  the  isochoric  process.”  But  it  is  not  always  so  obvious.  For 
instance,  the  quantity  L  =  T  (  // ,  —  r/2)  is  obviously  a  state  function — it  is  defined  only  in  terms  of 
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state  variables,  without  any  mention  of  the  particular  process  of  phase  transformation.  However, 
the  very  name  “the  latent  heat  of  phase  transformation”  says  to  us  that  it  is  the  rate  of  heat  supply 
in  a  certain  thermodynamic  process.  Unfortunately,  the  classical  textbooks  rarely  specify  clearly 
the  process  for  which  L  =  T  (rp  —  p2 )  is  the  rate  of  heat  supply.  The  lacking  analysis  creates  the 
impression  that  L  =  T  (rp  —  //2 )  describes  the  latent  heat  for  any  thermodynamic  process.  This  is 
what  practitioners  tacitly  assume  and  they  are  surprised  when  somebody  claims  that  this  is  a 
stretched  interpretation.  That  is  why,  in  this  report,  we  come  back  to  the  discussion  of  the  latent 
heat  concept.  We  deliberately  limit  ourselves  to  the  simplest  classical  one-component  system 
“liquid/vapor,”  which  everyone  has  learned  in  high  school  and  college. 

We  believe  that  the  general  situation  will  become  much  more  transparent  if  we  introduce  the 
concept  of  “the  latent  heat  capacity  L(t )  in  thermodynamic  process”  similar  to  the  concept  of  “the 
specific  heat  capacity  in  thermodynamic  process.”  Contrary  to  the  latent  heat  L  =  T  (rp  —  p2)  , 
the  latent  heat  L(t )  is  not  a  state  function.  However,  for  more  narrow  classes  of  thermodynamic 
processes,  it  leads  to  novel  state  functions.  In  addition,  we  introduce  a  novel  thermodynamic 
concept  of  “the  latent  work  in  thermodynamic  process,”  which  is  as  meaningful  as  the  concept  of 
“the  latent  heat  capacity  in  thermodynamic  process.” 


2.  Latent  Heat  and  Work  of  Thermodynamic  Process 


Consider  a  single-component,  two-phase  heterogeneous  system  “liquid/vapor”  (see  Fig.  1)  in  a 
pressurized  vessel.  If  the  piston  position  is  fixed,  then  the  phase  transformation  occurs  in  the 
isochoric  regime.  Also,  we  assume  the  system  can  be  thermally  isolated  (adiabatic  phase 
transformation)  or  placed  within  a  thermostat  (isothermic  regime.)  In  what  follows,  we  mark  as  1 
and  2  the  parameters  of  the  liquid  and  vapor  phases,  respectively. 

Consider  any  process  “x”  in  the  two-phase  heterogeneous  system  under  study.  Then,  depending 
on  the  pressure  and  thermal  loading,  the  external  heat  supply  Q,  pressure  P,  and  total  volume  Vtot 
become  functions  of  time.  The  same  is  true  regarding  the  thermodynamic  parameters  of  the 
phases  pi,  pi,  771,  ei,  ipi,  etc.,  and  p2,  p2,  r/2 ,  e2, ^2,  etc.  Generally  speaking,  thermodynamic 
processes  are  accompanied  by  the  evolution  of  the  total  mass  components  Mi,  M2. 

The  total  mass  of  the  substance  Mtot  is  assumed  fixed, 

Mi  (t)  +  M2  (t)  =  Mtot  ,  (6) 
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I  I  I 


Fig.  1  Fleterogeneous  system 

liquid/vapor  under  pressure 


or  in  differential  form 

dM\(t)  dM2(t ) 

dt  dt 


We  define  the  latent  heat  of  condensation  Lx  (t)  in  the  process  “x”  by  the  ratio 


Lx{t) 


dQ{t)/dt 

dMi(t)/dt 


(7) 


(8) 


Likewise,  we  define  the  latent  work  of  condensation  Wx  (t  )  in  the  process  “x”  by  the  ratio 


Wx(t)  =  P{t ) 


dVtot(t)/dt 
dM1  (t)  /  dt 


(9) 
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3.  Equations  of  Phase  Equilibrium 


Regardless  of  the  external  mechanical  and  thermal  loading,  the  following  phase  equilibrium 
conditions  should  be  satisfied: 

i)  the  condition  of  mechanical  equilibrium, 

Pi=P2  or  (p2ep)i  =  (p2ep)2  ;  (10) 

ii)  the  condition  of  thermal  equilibrium, 

Tx  =  X2  or  evi  =  en2  ;  and  (11) 

iii)  the  condition  of  equilibrium  with  respect  to  mass  exchange  between  phases, 

Xi  =  X2  or  (e  +  pep-r}ev)1  =  (e  +  pep  —  r}eTI)2  ,  (12) 

which  is  alternately  expressible  as 

{^  +  p/p)i  =  {^  +  p/p)2  ■  (13) 

The  quantity  x  is  known  as  the  chemical  potential  of  the  one-component  liquid  (or  gaseous) 
phase.  Therefore,  condition  12  can  be  called  the  chemical  equilibrium  condition,  although  no 
chemical  reactions  are  actually  occur  in  the  transformation  under  consideration. 

Also,  the  mass  conservation  condition  should  be  satisfied  for  any  external  loading.  All  additional 
conditions  depend  on  the  specifics  of  the  external  loading  associated  with  process  x. 


4.  Vaporization/Condensation  under  Fixed  Pressure  and  Prescribed  Heat 
Supply 


When  external  pressure  is  fixed  p  =  p0  the  vaporization/condensation  process  can  be  controlled 
by  the  external  heat  supply  Q(t). 

Differentiating  the  phase  equilibrium  Eqs.  10-12  along  the  process  trajectory  allows  us  to  obtain 
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governing  relationships.  First,  Eq.  10: 


(  2  dri  /  \  dP\  (2  drl  (  \  dP\  dPext 

e„-  +  p  („e)pp  -  j  i  =  (p  e„-  +  p  (pe)„  -  j  ^  =  — 

where  pext  is  the  constant  external  applied  pressure.  Equation  14  implies 

(  pepv  \  drj 2 


=  0 


(14) 


dpi 

dt 


pepri  \  drji 
(pe)„„/.  dt 


pp'  1 


dp2 

dt 


\Pe)pPJ  2  dt 


(15) 


Then,  differentiating  the  condition  of  thermal  equilibrium,  Eq.  11,  we  get 

drj  dp N 


^dt+eripdt 


drj  dp 
^dt+CrtpTt 


(16) 


Finally,  differentiating  the  condition  of  chemical  equilibrium,  Eq.  12,  we  get  after  some  algebraic 
manipulations 

(( eP)pP  -  Vevp)  %  +  (ewP  -  Pem)  ^ 


dt 


( (eP)pp  ~  VCvp)  %  +  ( eVPP  ~  r}evv) 


1 

dp 


(17) 


dt  y  "pr  '  dt 

The  system  given  by  Eqs.  15-17  can  be  rewritten  in  the  following  matrix  form: 


[G, 


dpi/dt 

O" 

drj  1  /  dt 

0 

dp-ifdt 

0 

drj2  /  dt 

0 

(18) 


where  Gp  is  the  following  4x4  matrix: 


[Gp]  = 


(p  ( Pe)  pp )  1 
0 


€"r]pl 


(p  epp) 1 

0 

Crjrj  1 


0 


(p  (pe 


'ppj  2 

6 r/p2 

(( eP)PP  -  VevP 


(( ep)pp  ~  pevp)  i  (evpp  -  r}evv)1 

The  determinant  of  the  matrix  Gp  is  given  by  the  following  formula: 


0 

(p  epv) 2 

^rjrj2 

( eVpP  —  Pew)s 


(19) 
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IGJ  = 


(pipe)  ) 

V  jppJ  1 

(p2ePri)i 

0 

0 

0 

0 

(pipe)  ) 
v  jppJ 2 

(p  epv) 2 

Hripl 

&qql 

dr]p2 

^?7772 

0 

0 

(jl^"qp)2  7l^-»7p2  (V^irq)  2  Vl^"t]r)2 

-pn) 


=  {vi  -  m)  (p(pe)wew  -  {pepnf)i  (p(pe)ppevv  -  (pep 

Introducing  the  specific  volume  v,  instead  of  p,  as  v  —  1/p,  we  can  rewrite  Eq.  20  as 

|Gp|  =  (r/i  '^2)  tq  (^Cvv^-rjr]  ^  ^2  (^vv&Tfr} 


,2 

^V1)J  2 


According  to  the  classical  thermodynamic  inequality,  we  have 


(20) 


(21) 


Gvv&rjrj  ^ 


Since  r/i  —  772  7^  0,  we  conclude  that 


|GP|^0 


Finally,  since  Gp  7^  0,  system  18  implies 


dpi  dpi  dp-2  dp  2 
dt  dt  dt  dt 


=  0 


(22) 


(23) 


According  to  Eq.  23,  because  two  independent  thermodynamic  parameters,  (pippi)  and  (p2,  P2), 
within  each  of  the  phases  remain  unchanged  in  the  processes  with  fixed  external  pressures,  all 
other  thermodynamic  parameters,  including  the  absolute  temperature,  remain  unchanged  as  well. 
This  fact  is  well  known  in  classical  thermodynamics. 

Let  us  now  calculate  the  latent  heat  Lp  for  the  isobaric  vaporization/condensation.  The  total 
internal  energy  E  and  the  total  volume  V  of  the  system  can  be  presented  as  follows: 


E  =  M1e(p1,p1)  +  M2e(p1,pi)  =  MIe(pI,pI)  , 

1= 1,2 


Vtot(t)  =  Vi(t)  +  V2(t)  = 


Mi(t)  M2(t)  x.  Mj 


+ 


1  pi 


Pl  P2 

Note  that,  in  light  of  Eq.  6,  Eq.  25  may  be  alternately  expressed  for  this  two-phase, 


(24) 

(25) 
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single-component  system  as 


VW t)  =  M1(t)  (---)+  — 
\P  1  P‘2  P2 


(26) 


Since  all  the  thermodynamic  parameters  of  the  phase  remain  unchanged  over  time,  we  can  employ 
the  First  Law  of  Thermodynamics  to  obtain  the  heat-transfer  rate  for  this  process  as  follows: 


dQ 

dt 


as  it  is  implied  by  the  following  chain: 

dQ  _  dE 
dt  dt 


( 1 

Ml 

ei  -  e2  +  p0  - 

-  - 

L  vpi 

P2/J 

dM\ 

dt 


(27) 


Po- 


dVt 


tot 


dt 


d  ^  „  r  .  ,  d,  ^  M/ 

-YiM,e{p,,V,)+po-R'£  — 


x  _  dMj  1  dMj 

,  .  dMx  (  1  1  \  dMi 

ei-e2  —  +p0 - —vr 

at  \p  i  p2  J  dt 


1  1  N 

ei  -  e2  +  po\ - 

\Pl  P2, 


where  we  have  relied  upon  relationships  6  and  24-26. 


dMx 

dt 


Comparing  the  general  definition  of  the  latent  heat  of  thermodynamic  process  with  relation  27  for 
isobaric  process,  we  arrive  at  the  following  relationship  of  Lp 


,  _  dQ  _  ,  /  1  1' 

Lp  —  -jvy-  —  ei  —  e2  +  p0 - 

dM i  VPi  p2, 


(28) 


Using  the  chemical  condition  of  phase  equilibrium,  Eq.  12,  we  can  rewrite  in  the  following 
classical  form: 

Lp  =  T0  (pi  -  r/2)  .  (29) 


For  the  latent  work  of  vaporization/condensation  at  fixed  external  pressure  we  get,  according 
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to  definition  9 


Wp(t)  =  p0 


dVtot(i)/dt 

dMi(t)/dt 


=  Po 


(30) 


where  the  time -derivative  of  total  volume  has  been  obtained  by  taking  the  derivative  of  Eq.  26.  In 
light  of  Eq.  13,  the  formula  of  the  latent  work,  Eq.  30,  can  be  rewritten  in  terms  of  the  free  energy 


jump: 


Wp(t)  =  po 


dVtot(t)/dt 

dMi(t)/dt 


fp2  -  i>l 


(31) 


Now,  we  can  rewrite  the  classical  Clausius-Clapeyron  formula,  Eq.  3,  in  the  following  transparent 
and  elegant  form: 


dp(T )  _  rji  -  r]2  _  Lp  po 

dT  ~  -P2X~  Wv  To 


(32) 


dlnp  _  To  Vi  -  rj2  _  Lp_ 
d  hi  T  po  pi1  -  P21  Wp 


5.  Vaporization/Condensation  at  Fixed  Total  Volume 


When  heating  under  fixed  volume,  the  phase  transformation  temperature  and  pressure  grow.  The 
latent  work  of  phase  transformation  Wv  vanishes.  The  formula  of  the  latent  heat  changes  also. 
Taking  the  time-derivative  of  Eq.  25,  we  obtain 

dV tot  _  dVj  _  (  1  dMr  _  Mj  dp  A  _  Q 

dt  j  dt  i\Pi  dt  pj  dt  J 


while  the  total-mass  M  conservation  conditions  implies 


E 


dMr 

dt 


0 


(35) 


For  the  heat  supply,  we  get  the  following  relationship: 


d,Q 

dt 


T  (771  -  V2) 


dMr 

dt 


E  MiT 


dpi 

dt 


(36) 
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as  it  is  implied  by  the  chain 


d,Q  _  dE 
dt  dt 


d  ^  ,  r  .  .  dM\ 

St?  M/e,  =  (e,  -  e2)  — 


V-  A/r  (  dPl  , 


=  (ei  -  e2) 
=  (ei  -  e2) 


rfMl  j_  i\/r  (  1  dPi  ,  rj.d'Pi' 

+  }^Mi  \  P  ■>  ..  +T- 


dt 

dM1 

dt 


pj  dt  '  dt 


+  £ 


p  dMj 


dpi' 


,  +  m,t 

Pj  dt  dt 


p  p\  dMi  _  dip 

=  I  ei  -  e2  +  —  -  —  )  -Jr-  +  2^  -^r 


Pi  P2  /  dt 


The  phase  equilibrium  conditions  imply  three  conditions: 

i)  the  mechanical  equilibrium  condition 


dp 

dt 


2  dr/ \  //  2  \  dp  2  dp\ 


2  \  dp  ,  . 

pep)p-  +  pew-]  = 


(37) 


ii)  the  thermal  equilibrium  condition 


dT 

dt 


'e  —  +  e  *L\  =  (e  ^  +  e  ^  •  and 

v  w  dt  +  w  dt  J  ,  ^  w  dt  +  w  dt  J  2  ’ 


(38) 


iii)  the  “chemical”  equilibrium  condition 


dx 

dt 


(( Pe)pp  ~  Vevp)  ^  +  (pevp  -  pevv) 


dt 


(( Pe)pp  ~  VeVP)  ^  +  (peVp  -  pem) 


dt 


dp 

dt 

dp 

dt 


J  i 


J  2 


(39) 


The  heat  flux  relation,  Eq.  36,  can  be  rewritten  as 


.  .  dMl  /i,r  dVl  1  dQ 


(40) 
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The  variable  dM2/dt  can  be  eliminated  between  Eqs.  34  and  35  to  read 

MT  dpi  +  (l_  _  l\  dM\  _  Q 
j  P2i  dt  \p2  pi)  dt 


In  matrix  form,  the  system  of  equations,  Eqs.  37-41,  may  be  expressed  as 


dpi/ dt 

0 

dip  /  dt 

0 

dp2/dt 

= 

0 

dp  2  /  dt 

(1/T)  •  (dQ/dt) 

dMi/dt 

0 

where  the  matrix  [Gyl  is  defined  as 


P  (Pe)PP 1 

P2epv 1 

~P  (Pe)PP 2 

- P2epri2 

0 

&r]pl 

67777 1 

&r)p2 

677772 

0 

')  PP~  '^ew)  1 

( Pevp  ~  Vevv)i 

-{(Pe)PP~Vevp)2 

—  ( Pevp  ~~  Vevv)2 

0 

0 

Mi 

0 

m2 

Vi 

— 

72 

Ml/pi 

0 

M2/p\ 

0 

P2 

— 

J_ 

pi- 

Thus,  with  Cramer’s  Rule,  we  arrive  at  the  following  formula: 


dMx  _  |Uy |  1  dQ 

dt  I  Gy  I  T  dt  ’ 


where  the  matrix  [Uy]  is  defined  as 

P  ( Pe)PPi  P2ePv i 

^"npl  ^riril 

[Uy]  =  .  x 

( (Pe)pp  -  Vevp)  1  ( Pevp  -  vemh 

My/p\  0 


-P  ( Pe)pp2 

^ rjp2 

’(( Pe)pp-Vevp 
M2/pl 


-P2epv  2 


^r/772 

{P^pp  ^^7777)^ 

0 
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Thus,  we  arrive  at  the  following  formula  of  the  latent  heat  at  fixed  volume: 


Lv  =  -T 


6.  Adiabatic  Vaporization/Condensation 


In  the  case  of  adiabatic  vaporization/condensation,  the  latent  heat  is,  by  definition,  equal  to 
Per  Eq.  36,  the  adiabaticity  condition  reads  (i.e.,  when  -'/'1 '  -u  ~  m 


dQ/dtt  =  0) 


dM\  x  '  dr] i 

=0-  - 


In  order  to  alternately  derive  Eq.  47,  we  may  first  obtain 


dMl  dVu 

— -p— 


]Ml  dX tot  rp  l\,r  d'dl 

— - p — - - hT)  Mi— - 

dt  dt  r  dt 


as  implied  by  the  chain 

dE  d  x  -  /  x  \  '  ( dMT  de(pi,rji) 

M=dt^  Ml  e(p •"  Pl)  =  ?  +  Ml  ~a — 


i  ^  1  ^  dMi  \ '  dVj  \ '  drji 

?i-^+P(---)\^r-pZir  +  TZM,w 


rp,  x dMl  dVt, 

where  Eq.  34  has  notably  been  used. 


dt  dt  r  dt 


Now,  combining  the  above  with  the  First  Law  of  Thermodynamics,  we  get 


dQ  dE  dX tot  rp  /  \  dMi  x.  -x  dp  j 

0  =  s  =  a + v^r = T  <"■  -  ^  nr+T^Mlit 
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which,  when  divided  by  T,  yields  Eq.  47. 


For  the  latent  adiabatic  work,  we  get  the  formula 


dAn 

dt 


\  '  Mi  dp i  \  '  ,  ,  dr) i  d ,\l < 


pj  dt 


(49) 


by  noting  that  dAv/dt  =  —dE/dt,  as  given  in  the  third  line  of  the  prior  chain. 


The  phase  equilibrium  conditions  still  imply  the  following: 

i)  the  mechanical  equilibrium  condition 


dp 

dt 


((p2ep\dt +p2epvdt)1~  (Vep) 


dp  2  drf 
p  dt  P  pn  dt  , 


(50) 


ii)  the  thermal  equilibrium  condition 


dT 

dt 


( e  —  +  e  ^  =  (e  ^  +  e  •  and 

K  vp  dt  +  m  dt )  ,  \  np  dt  +  m  dt)2  ’ 


(51) 


iii)  the  “chemical”  equilibrium  condition 


dx 

dt 


(( Pe)PP  ~  Vevp)  37  +  ( PevP  ~  Vevv)  37 


dr) 


dt 


(( Pe)PP  ~  Vevp)  %  +  ( PevP  ~  Vevv)  37 


dt 


dr) 


dt 


dt 


(52) 


The  system  of  equations  given  by  Eqs.  47  and  49-52  can  be  rewritten  in  the  following  matrix 
form: 


dpi/dt 

0 

dpi  /  dt 

0 

dpi/dt 

= 

0 

d/rp/dt 

0 

dMi/dt 

— dA^/dt 

(53) 
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where  [G^]  is  the  matrix 


[Gjj]  = 


p  ( Pe)ppl 

P2eprii 

~P  (Pe)PP 2 

-p2epr)  2 

0 

£77771 

6 )yp2 

^77772 

0 

')PP~VeVP)1 

(Per)p  ~  1leiTr])i 

-((Pe)pp  -  Vevp)2 

—  {Pevp  —  Vevn)  2 

0 

0 

Mi 

0 

m2 

C3 

1 

to 

pM1/p\ 

M\T 

pM2/pl 

m2t 

ei  -  e2_ 

(54) 


Thus,  we  arrive  at  the  following  formula: 


dM\ 

dAv  U,? 

dt 

dt  |  Gq  | 

where  U,?  is  the  matrix 

p  ( Pe)ppl 

P2epp 1 

-P  ( Pe)pp2 

-p2epv  2 

[u„]  = 

^77pl 

67^77! 

Cpp2 

^77772 

(( Pe)pp~Verip)1 

{per]p  ~  Vevv)i 

-{(Pe)pp-Vevp)2 

~{per/p  —  Verrn)  2 

0 

Mi 

0 

m2 

(55) 


(56) 


Thus,  the  specific  latent  work  Wv  of  vaporization/condensation  for  an  adiabatic  process  is  given 
by  the  formula: 


Wv  = 


|Gr 


IU, 


(57) 


7.  Conclusion 


We  introduce  several  novel  notions  relating  to  the  first-order  phase  transformation  in  liquid/vapor 
one-component  systems.  In  addition  to  the  standard  state  function  L,  called  the  latent  heat,  we 
introduce  the  concept  of  the  latent-heat-for-the-process  function  Lx.  Contrary  to  the  /.-function, 
the  /^-function  is  not  a  state  function  but  a  function  that  changes  depending  on  the  particular 
process  under  study.  The  function  Lx  coincides  with  L  only  for  the  phase  transformation  under 
fixed  pressure  and  temperature. 
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We  introduce  the  concept  of  the  latent  work  W  of  the  phase  transformation.  Here  again,  we 
distinguish  between  the  latent  work  W  as  the  state  function  and  the  latent-work-for-the -process 
function  Wx. 

Importantly,  this  report  lays  out  a  generalized  methodology  for  obtaining  the  latent  heat  and  latent 
work  for  variously  constrained  thermodynamic  processes.  Our  methodology,  though  lengthier 
than  the  textbook  derivation,  is  technically  straightforward  and  can  be  used  in  different 
circumstances. 

We  calculate  the  Lx  and  Wx  function  for  variously  constrained  thermodynamic  processes 
accompanied  by  phase  transformations.  For  the  classical  vaporization  process  under  fixed 
pressure  and  temperature,  we  give  a  new  elegant  form,  Eq.  33,  of  the  celebrated 
Clausius-Clapeyron  relationship. 
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Appendix.  A  Curious  Relation 


Curiously,  Holman3  derives  the  relationship  for  a  one-component  single-phase  liquid  as 


dp  drj 
dT  dv 


(A-l) 


which,  to  some  extent,  is  reminiscent  of  the  Clausius-Clapeyron  equation.  He  begins  with  the 
Helmholtz  free  energy  w  =  e  —  Tq,  from  which  an  exact  differential  may  be  taken: 


dip  =  de  —  T  drj  —  rj  dT 


(A-2) 


The  increment  of  internal  energy  de  may  be  eliminated  by  way  of  the  heat/work  energy  balance, 
de  =  T  dq  —  pdv,  where  v  is  the  specific  volume,  equal  to  1/p.  Substitution  of  this,  in  turn, 
eliminates  the  term  T  dq,  leaving  the  expression 


dip  =  —  p  dv  —  q  dT 


(A- 3) 


The  implication  of  this  relationship  is  that 

dip(T,v)  dip(T,v) 

- p= 

Furthermore,  one  may  equate  the  mixed  partial  derivatives  as  follows: 

d2ip  dp(T,v )  dq(T,v) 

dT  dv  dT  dv 


(A-4) 


(A- 5) 


The  similarity  of  Eq.  A- 5  to  the  Clausius-Clapeyron  relation  of  Eq.  3  is  obvious,  upon  realizing 
that  (r/i  —  p2)/(pr1  —  p^1)  is  simply  Aq/Av. 
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Intentionally  left  blank. 
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